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Abstract 

Surface mapping plays an important role in geometric processing. They induce both area and angular distortions. 
^^^^ If the angular distortion is bounded, the mapping is called a quasi- conformal map. Many surface maps in our 

^^ physical world are quasi-conformal. The angular distortion of a quasi-conformal map can be represented by Beltrami 

^^ differentials. According to quasi-conformal Teichmiiller theory, there is an 1-1 correspondence between the set of 

^S| Beltrami differentials and the set of quasi-conformal surface maps. Therefore, every quasi-conformal surface map 

Q can be fully determined by the Beltrami differential and can be reconstructed by solving the so-called Beltrami 

OJ equation. 

^\ In this work, we propose an effective method to solve the Beltrami equation on general Riemann surfaces. 

^^ The solution is a quasi-conformal map associated with the prescribed Beltrami differential. We firstly formulate a 

discrete analog of quasi-conformal maps on triangular meshes. Then, we propose an algorithm to compute discrete 
I— I quasi-conformal maps. The main strategy is to define a discrete auxiliary metric of the source surface, such that 

"^ the original quasi-conformal map becomes conformal under the newly defined discrete metric. The associated map 

^^ can then be obtained by using the discrete Yamabe flow method. Numerically, the discrete quasi-conformal map 

converges to the continuous real solution as the mesh size approaches to 0. We tested our algorithm on surfaces 

scanned from real life with different topologies. Experimental results demonstrate the generality and accuracy of 

our auxiliary metric method. 
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I. Introduction 
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QQ Teichmller theory 

^^ 

Xf^ Mapping between surfaces plays a fundamental role in digital geometric processing. In general, surface 

O mappings introduce distortions, which can be classified as area distortion and angular distortion. Map- 

^ pings without angular distortions are called conformal mappings. Last several years, there has been fast 

I* development of various techniques for computing conformal mappings and their applications in geometric 

K~A pi'ocessing. However, conformal mappings are not common in practice. Many mappings in our physical 

r\ world are quasi-conformal, which introduce bounded angular distortion. For example, deformations of 

c3 elastic shapes are quasi-conformal, such as human expression change, deformations of human organs, etc. 

In order to model surface mappings in the real world more effectively, it is crucial to study quasi-conformal 

mappings which allow for a much wider domain of applications. 

The theory of quasi-conformal mappings is nearly 80 years old and has been firstly studied by Ahlfors 
f\\, Grotzch [2J, Morrey [3 J and Lavrentjev [4J. Quasi-conformal mappings can be viewed as a gen- 
eralization of conformal mappings. Figure [T] illustrates the difference between a conformal map and a 
quasi-conformal map. Angular distortion can be characterized in the following intuitive way. Geometrically, 
a conformal mapping maps an infinitesimal circle on the source surface to a circle on the target surface, 
as shown in the first row. A quasi-conformal mapping maps an infinitesimal circle to an ellipse, as shown 
in frame (E) and (F). The eccentricity and the orientation of the ellipse can be represented by a complex 
valued function, the so-called Beltrami coefficient //. Specifically, the ratio between the two axes of the 
ellipse is given by ^^rr| , and the orientation of the axis is related to argfi{z) (see figure 2). 

W. Zeng and L.M. Lui contributed equally for this work 
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Fig. 1 

CONFORMAL AND QUASI-CONFORMAL MAPS FOR A TOPOLOGICAL DISK. (A) IS THE ORIGINAL FACE. (B) IS THE CONFORMAL 

MAPPING OF (A). (C) IS THE CIRCLE PACKING INDUCED BY (B). (D) IS THE CHECKERBOARD TEXTURE INDUCED BY (B). (E) IS THE 

QUASI-CONFORMAL MAPPING. (F) IS THE CIRCLE PACKING INDUCED BY (E). 



Beltrami coefficient is defined on a local chart. Globally, Beltrami coefficient is represented by the 
Beltrami differential, which is independent of the choice of local parameters. According to quasi-conformal 
Teichmiiller theory, for general surfaces, there is an one to one correspondence between the set of quasi- 
conformal maps and the set of Beltrami differentials. In other words, every quasi-conformal map can be 
fully determined by the Beltrami differentials and is unique up to a finite dimensional group. Conversely, 
given a particular Beltrami differential //(z)^, we can reconstruct the quasi-conformal maps associated 
to /x(z)^. The Beltrami differential captures the most essential information of the surface mappings. 
Therefore, by adjusting /i(z)^, we can reconstruct a surface mapping with desired properties. 

Quasi-conformal mappings have been studied extensively in complex analysis [IJ, [2J, [5], ^6}. Ap- 
plications can be found in different areas such as differential equations, topology, Riemann mappings, 
complex dynamics, grid generation and so on Q, [H, M, M, O, [[13, IE]], CHl. Despite the 
rapid development in the theory of quasi-conformal geometry, the progress on computing quasi-conformal 
mappings numerically has been very slow. In fact, developing an effective numerical algorithm to compute 
quasi-conformal mapping remains a challenging problem. 

Recently, there has been a few work on numerical quasi-conformal mapping techniques on the complex 
plane based on solving differential equations with finite difference or finite element methods. Most of 
these methods deal with simple domains in the complex plane and cannot be applied on arbitrary regions. 
Furthermore, to the best of our knowledge, no work has been done on solving the Betrami equation on 
general Riemann surfaces. In this work, we are interested in developing an effective numerical algorithm 
to compute the quasi-conformal mapping on general Riemann surfaces of any genus. Of course, the 
developed algorithm could be easily applied to any arbitrary regions in the complex plane C, since any 
regions in C are Riemann surfaces. 

The fundamental problem in this paper is to find the quasi-conformal map 0, associated to a given 
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Fig. 2 
Illustration of how the Beltrami coefficient /h measures the distortion by a quasi-conformal mapping that maps a 

small circle to an ellipse with dilation k. 



Beltrami differential. This can be done by solving the Beltrami equation using the proposed auxiliary 
metric method. We firstly formulate a discrete analog of the quasi-conformal map on triangular meshes. 
Then, we propose an algorithm to compute the discrete quasi-conformal map. The basic idea is to construct 
a discrete auxiliary metric based on the given Beltrami differential, such that (j) becomes a conformal map 
under the auxiliary metric. We then use the Yamabe flow method to compute the conformal map under 
the new discrete metric. The resulting map is the desired quasi-conformal map. Numerically, the discrete 
quasi-conformal map converges to the continuous real solution as mesh size tends to 0. 

The paper is laid out in the following way: Section II briefly review the most related works in the 
field; Section III introduces the theoretical background; Section IV describes the discrete Euclidean and 
hyperbolic Yamabe flow method; Section V explains how the Beltrami equations can be solved on general 
Riemann surfaces using the auxiliary metric. Section VI focuses on the computational methodologies; 
Section VII reports the experimental results; the paper is concluded in Section VIII. 

II. Previous work 

Conformal mapping has been broadly applied to surface parameterizations in digital geometry process- 
ing. Here we only review the most related works. We refer readers to the thorough surveys of ifTSll . lfT6l . 
^n\ for various kinds of mesh parameterization techniques. 

Levy et al. in ifTSll applied the Cauchy-Riemann equation for mesh parameterization and provided 
successful results on the constrained 2D parameterizations with free boundaries. Desbrun et al. in [ ?9l 
minimized the Dirichlet energy defined on triangular meshes for computing conformal parameterization. 
Angle based flattening method (ABF) was introduced in [l20l|, fl\\. Linearized version of ABF is introduced 
in II22I . Gu and Yau in ll23]| computed the conformal structure using Hodge theory. Gortler et al. in [l24l 
used discrete 1 -forms for mesh parameterization with several holes. Ray et al. in ||25]| used the holomorphic 
1-form to follow the principle curvatures for the quad remeshing purpose. Kalberer et al. f26l use branched 
covering to convert a given frame field on the surface to a vector field on the covering space. Spherical 
parameterizations are introduced in ll27ll . [|28]| . High genus surface parameterization is pioneered by Grimm 
and Hughes in ll29ll . Recently, hyperbolic parameterization is introduced in ll30ll . 

Circle pattern was proposed by Bowers and Hurdal [31], and has been proven to be a minimizer of 
a convex energy by Bobenko and Springborn [32J. An efficient circle pattern algorithm was developed 
by Kharevych et al. [331 . Discrete Ricci flow was introduced by Chow and Luo in [l34l and applied to 
graphics in ll30ll . Ben-Chen et al. introduce an efficient method for scaling metrics to prescribed curvatures 
in Hal. 
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The theory for combinatorial Euchdean Yamabe flow was introduced by Luo in ll36l . The theory for 
hyperbohc curvature flow was introduced in 071 . Springborn et al. Il38]| identifies the Yamabe energy with 
the Milnor-Lobachevsky function and the heat equation for the curvature evolution with the cotangent 
Laplace equation. 

Recently, there has been various work on numerical quasi-conformal mapping techniques based on 
solving elliptic equations in the real plane with finite difference or finite element methods. Using finite 
difference methods to compute quasi-conformal maps on complex plane were proposed by Belinskii et al. 
[[9J and Mastin and Thompson [14J. These methods are difficult to implement for arbitrary regions. A finite 
difference scheme for constructing quasi-conformal mappings for arbitrary simply and doubly-connected 
region of the plane onto a rectangle was developed by Mastin and Thompson ll39l . Vlasynk HOll applied 
similiar techniques for mappings of doubly connected and triply connected domains onto a parametric 
rectangle. A finite element based method was implemented by Weisel [41 J. In [42J Daripa proposed a 
numerical construction of quasi-conformal mappings in the plane using the Beltrami equation. The author 
presented an algorithm for the evaluation of one of the singular operators that arise in solving the Beltrami 
equation. The author subsequently applied the same method for numerical quasi-conformal mappings of 
exterior of simply connected domains onto the interior of a unit disk using the Beltrami equation [13]. This 
method was further extended to the quasi-conformal mapping of an arbitrary doubly connected domain 
with smooth boundaries onto an annulus Qr = {a : R < a < 1} [43J. All of these methods deal with 
simple domains in the complex plane and cannot be applied on arbitrary regions. Furthermore, as far as 
we know, no work has been done on solving the Betrami equation on general Riemann surfaces. In this 
work, we are interested in developing an effective numerical algorithm to compute the quasi-conformal 
mapping on general Riemann surfaces of any genus. 

III. Theoretical Background 

In this section, we briefly introduce the major concepts in differential geometry and Riemann surface 
theory, which are necessary to explain the quasi-conformal maps. We refer readers to ll44l , ll45il for detailed 
information. 

A. Beltrami Equations and quasi-conformal maps 

Let / : C ^ C be a complex function. The following differential operators are more convenient for 

d Id . d d Id . d 
dz ' 2 dx dy ' dz ' 2 dx dy ' 

f is said to be quasi-conformal associated to /x if it is orientation-preserving and satisfies the following 
Beltrami equation: 

where /x(z) is some complex- valued Lebesgue measurable function satisfying sup |//| < 1. /x is called the 
Beltrami coefficient of /. The Beltrami coefficient /x gives us all the information about the conformaity 
of/ (See Figure [2]). 

If ii{z) = everywhere, / is called holomorphic. A holomorphic function satisfies the well-known 

Cauchy-Riemann equation 

df 

dz 

Suppose S' is a surface embedded in M^, with the induced Euclidean metric g. Let [/^ C S' be an open 
set on S, with local parameterization is (l)a • Ua ^ C, such that the metric has local representation 

g = e^^^'^Uzdz, 
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Fig. 3 
The figure illustrates the definition of quasi-conformal maps between Riemann surfaces 



then {lJo,,(\)a) is called an isothermal coordinate chart. We can cover the whole surface by a collection of 
isothermal coordinate charts. All isothermal coordinate charts form a conformal structure of the surface. 
The surface with a conformal structure is called a Riemann surface. 

Suppose Si and 5^2 are two Riemann surfaces. ([7^, (j)a) is a local chart of Si, (V^, 7/^/3) is a local chart 
of 5^2. / : S'l ^ 5^2 is a conformal map if and only if 

is bi-holomorphic for all 0^ and iIj^. For simplicity, we still use / to denote its local representation. Then 
a conformal map / satisfies |f = 0. 

The definition of quasi-conformal maps of plane domains can be extended to Riemann surfaces. Instead 
of using the Beltrami coefficient, a global quantity called Beltrami differential is used, which is independent 
of the choice of local parameters. 

Definition 3.1 (Beltrami Differentials): A Beltrami differential m(^)^ on a Riemann surface R is an 
assignment to each chart z^ on Ua an Loo complex- valued function /x"^ defined on Za{Ua) such that: 



/.«(..) = /M(^)/(^) 



KAj/Cify KAJ/Cif^ 



(2) 



on the region which is also covered by another chart Zj^. 



Now, a quasi-conformal map between Riemann surfaces can be defined as follow (see Figure |3] for the 
geometric illustration): 

Definition 3.2 (quasi-conformal maps between Riemann surfaces)^ An orientation-preserving home- 
omorphism / : aSi ^ 6^2 is called quasi-conformal associated to l^{z)^ if for any chart Za on Si and any 
chart Wis on S2, the map /^^ \= Wjs o f o z~^ is quasi-conformal associated with ij!^{za) 

Note that the above definition is well defined. On a region of Si covered by two different charts z^ 
and Zc,/, we have 



M 






dz'^ dz'^ dZa dZa' dZa dZa'' 



' dZfy,/ dZa' ' 
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This is guaranteed by Equation |2| Also, the definition does not depend on the chart w^ used in the range 
of /. Let Wis and Wis^ be two different charts on the range of /. We have 



df^^' df^^' _ dwp^df^^ dw^df^ dw^df^ dw^df^ 



dZa dZa 



dw^f 



since w^f is holomorphic and so ^^ = 0. 

With the above definitions, we can now formulate our problem of interest as follow: 

PROBLEM : Let Si and 5^2 be two Riemann surfaces with the same topology. Given a Beltrami 
differential /i(z)^, we are interested in finding the map / such that for each chart z^ on Si and each 
chart Wis on 5^2, the map /^^ \= Wjs o f o z~^ satisfies the partial differential equation: 

-^=^ = M ^a ^ — 

OZa OZa 

B. Ricci Flow 

Let S' be a surface embedded in M^ with the induced Euclidean metric g. We say another Riemannian 
metric g is conformal to g, if there is a scalar function u \ S ^^, such that g = e^^g. 
The Gaussian curvature induced by g is 

where Ag is the Laplace-Beltrami operator under the original metric g. The above equation is called 
the Yamabe equation. By solving the Yamabe equation, one can design a conformal metric e^^g by a 
prescribed curvature K. 

Yamabe equation can be solved using Ricci flow method. The Ricci flow deforms the metric g(t) 
according to the Gaussian curvature K{t) (induced by itself), where t is the time parameter 

^ = 2(K-K(t))^.,(t). 

The uniformization theorem for surfaces says that any metric surface admits a conformal metric, which 
induces constant Gaussian curvature. The constant is one of{ — 1,0,+1}, determined by the topology of 
the surface. Such metric is called the uniformization metric. Ricci flow converges to the uniformization 
metric. Detailed proofs can be found in fi46iil and HTll . 

IV. Discrete Euclidean and Hyperbolic Curvature Flow 

In practice, most surfaces are approximated by simplicial complexes, namely triangular meshes. Suppose 
M is trianglar mesh, V, E, F are vertex, edge and face set respectively. We use vi to denote the ith vertex; 
[vi, Vj] the edge from vi to Vj\ [vi, Vj.Vk] as the face, where the vertices are sorted counter-clock- wisely. On 
triangular meshes, we can derive a discrete version of Ricci flow, called the discrete Yamabe flow, which 
is analogous to the curvature flow on smooth surfaces. In this section, we describe in detail the discrete 
Euclidean and hyperbolic curvature flow that converge to the Euclidean and hyperbolic uniformization 
metric respectively. 

On the discrete mesh, we can define the discrete metric, which is similar to the Riemannian metric. 
Basically, the discrete metric gives the length of each edge. 
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Definition 4.1 (Discrete Metric): A discrete metric on a mesh M is a function I \ E ^ M+, such 
that on each triangle [vi,Vj,Vk], the triangle inequality holds, 

Discrete metric represents a configuration of edge lengths. As shown in figure |4]^ A), different background 
geometries can be assigned to a mesh. 

Definition 4.2 (Background Geometry): Suppose M is a mesh with a discrete metric. If all faces 
are Spherical (or Euclidean, or Hyperbolic triangle), then the mesh is with Spherical (or Euclidean, or 
Hyperbolic) background geometry, denoted as §^ (or E^, or H^). 

Discrete metric determines the comer angles on each face by the cosine law, 

r cos-i '^^ E2 

9- = < ^^^^^ (3) 

^ I ^^^.-1 cosh Zj cosh ^fc -cosh /^ ^2 ^ ^ 

l^ 2sinh/j sinh/fc 

The discrete Gaussian curvature is defined as angle deficient. 

Definition 4.3 (Discrete Gaussian Curvature): Suppose M is a mesh with a discrete metric, which 
either in Euclidean or hyperbolic background geometry, [vi, Vj, Vk] is a face in M, Oj^ represent the corner 
angle at Vi on the face. The discrete Gaussian curvature of Vi is defined as 

The total Gaussian curvature is controlled by the topology of the mesh, 

Y,Ki + xY,Ajk = 27Tx{M), (5) 

i ijk 

where Aijk is the area of the face [vi, Vj.Vk], A is zero if M is with Euclidean background metric, —1 if 
M is with hyperbolic background metric. 

Definition 4.4 (Discrete conformal deformation): Let X be a triangulation mesh. Suppose / and L 
are two different discrete metrics on K. We say L is a discrete conformal deformation of / if there exists 
a function u :V ^R, where V is the set of all vertices of K, such that for all edges [vi, Vj] on Ki: 

Lmv^)J{vj)]) = I 2sinh-^e-(-sinh(Ml)e-(^^- M^ ^^^ 

where [vi.Vj] is an edge on K. QB)) 



Note that the definition of discrete hyperbolic conformality was due to Springborn et al. OTll . u : V ^R 
is called the discrete conformal factor. The discrete conformal factor is a function defined on every vertices 
of K. 

Definition 4.5 (Discrete Yamabe Flow): Let K^ denote the user defined curvature at Vi, the discrete 
Euclidean and Hyperbolic Yamabe flow has the same formula 

dUo — 

-^ = K.-K,. (7) 

The discrete Yamabe flow iteratively modifies the discrete conformal factor and hence define a sequence 
of conformal deformation. The initial metric is deformed into the uniformization metric. 
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Vk 




hj = l(If(vi), f(vj)]) ° „^ Lij = L([f(vi), f(Vj)]) 

(B) 



Fig. 4 
(A) SHOWS THE Euclidean and Hyperbolic triangle. (B) shows the discrete surface Yamabe flow 



In both Euclidean and hyperbolic cases, one can easily observe that under the discrete conformal 
deformation (see equation [13]), 



de, 09, 



duj dui 



(8) 



Therefore, the differential 1-form "^.Kidui is a closed 1-form. The following Yamabe energy is well 
defined. 




Fig. 5 
The figure illustrates the idea of discrete curvature flow. A planar domain with non-zero curvature on the 

BOUNDARY IS BEING DEFORMED ITERATIVELY INTO A FLAT METRIC. UNDER THE DEFORMED METRIC, CURVATURE IS ZERO 
EVERYWHERE IN THE INTERIOR AND IS CONSTANT ON THE BOUNDARY (I.E. THE BOUNDARY BECOMES A CIRCLE). 



Definition 4.6 (Discrete Yamabe Energy): Let u = (t^i,t^2, • ' ' ^^n) b^ the conformal factor vector, 
where n is the number of vertices, Uq = (0, 0, • • • , 0), then the discrete Euclidean and Hyperbolic Yamabe 
energy has the same form as 

^(u)= r5^(z,-x,)rf^,. (9) 

The discrete Yamabe flow is the negative gradient flow of the Yamabe energy. The Hessian matrix (hij) 
of the discrete Euclidean and Hyperbolic Yamabe energy can be computed explicitly. 

Theorem 4.7: Let [vi.Vj] be an edge, connecting two faces [vi,Vj,Vk] and [vj,Vi,vi], then Hessian 
matrix (hij) satisfies 

/,.^^ + ?^. (.0, 



Hj 



duj duj 
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where the summation goes through all faces surrounding vi, [vi,Vj,Vk]. Here, 





A{cj+l){ck+l] 


1 


de, 


r cot^fc 


E2 


duj 


S Ci+Cj-Ck-l 
I A(cfc+1) 


H2 



a^. j - cot ^j- cot 4 E2 

^ ^ "N 2ciCjCk-c'j-c1-\-CiCj-\-CiCk-Cj-Ck ^2 (1^) 

(13) 

where in H^ case, c/e is C08h{yk) and ^4 is sin(0/e) sinh(^^) sinh(7/j). 

Proof: See Appendix. ■ 

By carefully examining the positive definiteness of the Hessian, the convexity of the Yamabe energy 
can be obtained. 

Theorem 4.8: The discrete Euclidean Yamabe energy is locally convex on the space of ^^t^i = 0. 
The discrete hyperbolic Yamabe energy is convex. 

Proof: See Appendix. ■ 

The discrete Yamabe energy can be optimized using Newton's method directly. Given the mesh M, a 
conformal factor vector u is admissible, if the deformed metric satisfies the triangle inequality on each 
face. The space of all admissible conformal factors is not convex in either Euclidean or Hyperbolic case. 
In practice, the step length in Newton's method needs to be adjusted. Once triangle inequality doesn't 
hold on a face, edge swap needs to be performed. After finitely many such surgery operations on the 
triangulation mesh, there will be no singularity developed in the normalized discrete Yamabe flow. We 
can then prove that the discrete Yamabe flow converges exponentially fast to the discrete metric with 
constant curvature. Specifically, we have 

Theorem 4,9: If no singularity develops in the discrete Yamabe flow after finitely many of surgery 
operations, the solution converges exponentially fast to a discrete metric with constant curvature as time 
approaches infinity. In other words, 

\K,{t)-K,\<c,e-'^' 

for some constants ci and C2 and Ki{t) is the discrete curvature at vertex Vi at time t. 

Proof: See Appendix. ■ 

Figure [6] shows the exponential convergence of the discrete Yamabe flow method. The Yamabe energy 
at each iteration of two different real human faces are plotted on the left and right respectively. 

V. Solving Beltrami Equations by Auxiliary Metric 

In this section, we prove the main theorem of this paper, which allows us to define an auxiliary metric 
to solve the Beltrami equation on the Riemann surface. It also lets us define the discrete version of 
quasi-conformal maps between meshes. 

Theorem 5.1 (Auxiliary Metric associated to a Beltrami differential): Suppose (5'i,gi) and (5^2, g2) 
are two metric surfaces, : S'l ^ 5^2 is a quasi-conformal map, the Beltrami differential is /x^. Let z 
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Fig. 6 

The figure shows the exponential convergence of the discrete Yamabe flow method. The Yamabe energy at each 

iteration of two different real human faces are plotted on the left and right respectively. 



and w be the local isothermal coordinates of Si and 5^2 respectively, then gi is e^^^^^^dzdz and g2 is 
e^^'^^'^^dwdw. Define an auxiliary Riemannian metric on Si, 



gi = e 



2\i{z) 



\dz + /xrfz|^. 



(14) 



The auxiliary metric gi is vv^ell-defined and the map : (aSi, gi) ^ {S2, g2) is a conformal map. 



Proof: We first prove the auxiliary metric gi is w^ell-defined. Consider the region w^hich is covered by 
tvv^o different charts z^ and z^. Suppose the local representations of gi under z^ and z^ are e^^^^^^rfz 

and e^^^^^^ dz^ dz(^ respectively. 



Since ^ = 0, w^e have 

dzP 



Also, 



dz- = '^dz^ + Hdz^ = ^ 

dzl^ dzP dzP 



.dz'^, 



dz^ 



g2A4z)^^a^^ = e^^^^'^^ldz''^' = e'^^^'>\^\''\dz 



a\2 ^ ^2Aa(2:)|:;^|2|^^/3|2 ^ g2A^(2:) |^^/3|2 

dzi^ 



This gives, e^^/^^^) = e^^-^^^lgl^. 
Thus, 



dzP dzP 

_ g2A„(z«)|^|2| > fl , „a/^ /f^N^-^|2 

_ ' Uz^^ '"^^ +^ ^dzPUzP'"^' ' 

= e2A,(./')|^^/3^^/3^^^|2 

To see the map : (S'i,gi) ^ (5*2, §2) is a conformal map, let (f)*g2 denote the pull back metric. 



Under the pull back metric, the map : {Si, (f)*g2) -^ ('S'2,g2) is isometric. 

= ^^dz + ^^c^ 
— -j){dz + iJ,dz). 



d(j){z) = ^dz + ^dz, 
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/^; 



3(^4 



ij 



dz 



Fig. 7 
The figure illustrates the idea of discrete local chart and discrete Beltrami differential. 



Therefore, 



^g2 



.^MHz)) 



d(l){z) 
dz 



^'\dz + iJ.dz\ 



According to the definition of gi in Equation 

gi. Because cj) : (S'i,(/)*g2) -^ (5^2, g2) is isometnc, therefore 



14 



g^ ^ g2A2Wz))-2Ai(z)|M£)|2g^^ 



:(5i,gi)^(52,g2: 



z!)*g2 is conformal to 
is conformal. ■ 



This theorem tells us in order to solve the Beltrami equation on the Riemann surface, we simply need 
to define a new auxiliary metric associated with the prescribed Beltrami differential. We can then solve the 
Beltrami equation by computing a conformal map associated with the newly defined metric. Specifically, 
we need to compute a conformal map cj) : (S'i,gi) -^ (5^2, g2) (from Theorem 5.2). It can be done by 
flattening (5'i,gi) and (5^2, g2) onto the unit disk. Let (fi : (5'i,gi) -^ D and 9^2 : (5^2,^2) ^ D be the 
conformal parameterizations of Si and S2 respectively. The conformal map : (ASi,gi) -^ (5^2, g2) can 
be obtained by: (/) = (^2^ oLpi. To compute the conformal parameterizations Lpi (i = 1, 2), it suffices to find 
the flat metric of (S'i,gi) and (AS2,g2). We can apply the discrete Yamabe flow method to conformally 
deform the original metric on the triangular mesh to the uniformization metric. 

This motivates us to give a discrete definition of quasi-conformal maps between triangular meshes. 

Definition 5.2 (Discrete Local Charts): Let X be a triangular mesh. A mesh K^ is called a submesh 
of K if every vertices, edges and faces of K^ belong to K. A discrete local chart z^ : K^ ^ z^{Ka) C C 
is a discrete conformal map from K^ to a mesh z^{K^) embedded in the complex plane (See Figure [t]). 

Since the triangular meshes we consider are discrete approximations of smooth surfaces, we can assume 
that the triangular meshes are covered by a collection of discrete local charts. 

Definition 5.3 (Discrete Beltrami Differential): A discrete Beltrami differential jjl^ is an assignment 
to each discrete local chart z^ on K^, din L^ complex- valued function /i^ defined on z^{Kc,) such that: 

(15) 



( 



/X"(^,) +/X"(^,)x^"(^j) - ^"(^^) _ //(^^) +M^(^j)x^^(^j) - Z^{^i) 



Z^{Vj) - Z^{Vi) 



z^{v,)-z^{v,) 



where \vi, v^ is any edge in the region which is also covered by another chart z^ (See Figure [?]). 



By letting /x^ 



Equation 15 can be simplified as: 



^"^^'^+^"^"''^ 4 = /Aliil+^AM, rf^a = ^a(^^.) _ ^a(^.) and dz% = z^{v,) - zP{v,), 
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dz^- p dz 



13 



which is an analog of Equation |2} 






Definition 5.4 (Discrete quasi-conformal Map): Let /i^ be a given discrete Beltrami differential. A 
map (j) : (Ki, /) -^ {K2, L) between meshes^i and K2 is called discrete quasi-conformal if: with respect 
to a new metric I on Ki, the map (j) : {Ki, I) -^ {K2, L) is discrete conformal where: 

l{[vi,Vj\) := l{[vi,Vj] y ^'^Uz'a \ dzij = z{vj) — z{vi); fiij = ^'^^^ for any local isothermal coordinates 
zof Ki. 

I is called the discrete auxiliary metric associated with /x^. Note that the definition is well-defined. 
Suppose an edge [vi, Vj] is covered by both charts z^ and z^, we have 



l{[v^,Vj]) 






^(N, 


^i]) 1 + 




^(bi, 


vj]) 1 + 




/(T., 


., n <• 


+/igd4 



After the discrete auxiliary metric / is computed for each edge, we use the discrete Yamabe flow 
method introduced in Section IV to deform the metric into the uniformization metric and obtain the 
quasi-conformal map associated with the given Beltrami differential. 

VI. Computational Algorithms 

In this section, we describe in detail the algorithms for Euclidean and Hyperbolic Yamabe flow. We also 
describe the algorithm for computing the discrete auxiliary metric associated uniquely with the prescribed 
discrete Beltrami differential. 

A. Yamabe Flow 

Algorithm 1 describes the details for Euclidean and Hyperbolic Yamabe flow. 

Algorithm 1 : Euclidean & Hyperbolic Yamabe Flow 
Require : A triangular mesh M, the target curvature K. 



1 

2 
3 
4 
5 

6 
7 
8 
9 
10 
11 



Initialize Ui to zero for each vertex Vi, 
Repeat 

For each edge, compute the edge length using formula |6} 

For each face [vi, Vj.Vk] compute the corner angles 9i, 9j and 9k. 

For each face [vi,Vj,Vk], compute |^, ^, and ^ using formula 



For each face [vi,Vj,Vk\, compute |^, ^ and ^ using formula 



Construct the Hessian matrix H using formula [TO] and [TT 

Solve linear system H5u = K — K. 

Update discrete conformal factor u ^ u + Su. 

For each vertex Vi, compute the Gaussian curvature K^. 

until max^.^M \Ki — Ki\ < e 
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B. Embedding 

After computing the discrete metric of the mesh, we can embed the mesh onto M^ or H^. Here, we 
discuss the hyperbohc case. Euchdean case is very similar. 

In this work, we use Poincare disk to model the hyperbolic space H^, which is the unit disk |^| < 1 

on the complex plane with the metric 

, 9 Mzdz 
ds^ - 



{i-zzy 

The rigid motion is the Mobius transformation 

Z - Zq 



z ^ e 



1 - z^z 
where Q and zq are parameters. A hyperbolic circle (c,r) is also a Euclidean circle (C, K) with 

^ 2 - 2/x2 



c; 



1 — /x^cc 



.2 r<r^ CC - //- 



1 — /x^cc' 

where /x = tanh |. Given a two points p and g on H^, the unique geodesic through them is a circular arc 
joining them and is perpendicular to the unit circle. 

Basically, we can isometrically flatten triangle by triangle using the hyperbolic cosine law, as described 
in algorithm 2. . 

Algorithm 2 : Embed on H^ 



Require : A triangular mesh M, a set of fundamental group generators intersecting only at the base 
point p, using the algorithm in ll30ll 

1) Slice M along the base loops to form a fundamental domain M. 

2) Embed the first triangle [^0,^1,^2] ^ ^, 

r(^o) = (0,0),r(^i) =tanh— ,r(^2) =tanh— e<'. 

3) Put all the neighboring faces of the first face to a face queue. 

While : the face queue is not empty 

4) Pop the first face from the queue \vi, v^,V]^ 

5) Suppose Vi and v^ has been embedded, compute the intersection of two hyperbolic circles 

by converting them to Euclidean circles, 

6) r(^/e) is chosen the keep the orientation of the face upward. 

7) Put the neighboring faces, which haven't accessed, to the queue. 
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C. Auxiliary Metric 

The following algorithm 3 computes the discrete auxiliary metric based on the discrete Beltrami 
differential. 

Algorithm 3 : Auxiliary Metric 

Require : A triangular mesh M, a conformal parameterization z : V ^ C, discrete Beltrami differential 

fi:V ^C. 

1) For each edge [vi, Vj], compute the edge length lij using the Euclidean metric 
For all : Edge [vi.Vj] do 

2) dz ^ z{vj) — z{vi) 

3) /x^ Uii{vi) + fi{vj)) 

\dz-\-/j,dz\ 



4) A^ 
lij <- 
end for 



\dz\ 



VII. Experimental results 



We implement our algorithm using generic C++ on Windows platform. The linear systems are solved 
using conjugate gradient method. The experiments are carried out on a laptop with 2.0 GHZ CPU, 3.00 G 
RAM. The human face surfaces are captured using phase shifting structured light method. Computational 
time is reported in Table [1} 

A. Quasi-conformal Map for Genus Zero Surfaces 

Figure [8] shows the experimental results for quasi-conformal maps of the human face surface. The 
original face is shown in the top left corner. Four corner vertices are selected on the boundary, shown 
as Po,Pi,P2,PS' We set the target curvature to be | for those corner vertices, and zero for all vertices 
everywhere else. The Yamabe flow conformally maps the surface onto a planar rectangle. The corner 
vertices are mapped to the rectangle corners. We set the left lower corner to be the origin, the edges to be 
parallel to the axes, the width to be 1. Then the height h gives us the conformal module of the original 
face surface with four fixed corners (topological quadrilateral). This provides us the conformal parameter 
z of the surface. 

In the second and the third row, we set different Beltrami coefficients. The image of the quasi-conformal 
map is shown on the left, the circle packing texture mapping is shown on the right. The Beltrami coefficient 
is set to be // = ^^^ , with different values of zq for different cases. It is obvious that, the conformal 
module of the surface changes with different Beltrami coefficients. 

Figure |9] and [To] show the quasi-conformal maps for genus zero surfaces with multiply holes. In Figure 



[9| the human face surface is sliced open along the lip of the mouth which results in a doubly-connected 
open surface. Again, we set different Beltrami coefficients and compute the associated discrete auxiliary 
metrics. Using the discrete Yamabe flow, we conformally map the surface onto the annulus with respect 
to different auxiliary metrics. The target curvature is set to be zero in the interior and constant along the 
boundaries. The radius of the inner circles are different with different Beltrami coefficients, indicating 



a change in the conformal module. Figure 10 shows the similar results for the genus zero human face 



surface with three slices (topological disk with 3 holes). 

B. Composition of Quasi-conformal Maps 

In the following experiment, we test the accuracy of our algorithm by computing the composed quasi- 
conformal maps using different approaches, and comparing their difference. If our method is accurate, 
the difference between the results obtained from the two approaches should be small. 
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u-iz) = 






zo = 



Z[) = 0.5 + /0.5/i 





Z() = I + ih 



Fig. 8 
CONFORMAL AND QUASI-CONFORMAL MAPS FOR A FACE SURFACE. THE CONFORMAL PARAMETER DOMAIN IS A RECTANGLE WITH 



UNIT WIDTH AND h HEIGHT ON THE IST ROW. FOR ALL OTHER ROWS /i : 



2^1 + h' 



-WITH DIFFERENT Zq'S. 



Let f : S ^ Di ht a quasi-conformal map with Beltrami coefficient fif, g : Di ^ D2 with Beltrami 
coefficient fig. Then the composed map g o f : S ^ D2 should have the Beltrami coefficient 



l^g^f ~ ^ ^ -TT- 



l + fif{figof)r' 



(16) 



where r = ^. 

Jz I I 

As shown in figure [TTJ in our experiment, the original surface is a human face surface with four corner 
points (a topological quadrilateral), as shown in (a). We compute its conformal parameter domain, as 
shown in (b), (c) and (d). Then we set the Beltrami coefficient fif to be 0.15 + i0.15, and use our method 
to compute a quasi-conformal map f : S ^ Di. The mapping result of / is shown in (e) and (f). We set 
fig = 0.15 + i0.15, and compute the quasi-conformal map g : Di ^ D2, as shown in (g) and (h). 

We use the formula in Equation [16] to compute the Beltrami coefficient for the composed map figof = 
0.34 + z0.12. We then solve the Beltrami equation h^ = figofh^ to get a quasi-conformal map h : S ^ D2. 
In theory, h should coincide with go f. Our experimental result shows that h is consistent with g o f. By 
comparing the result in (g) and (i), we can see the results ofgof and h ait almost identical. We further 
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,«(z) = 0.2 




Fig. 9 

CONFORMAL AND QUASI-CONFORMAL MAPS FOR A TOPOLOGICAL ANNULUS. 



measure the deviation between them numerically, using the following formula, 

1 



d{f,g) 



f Ifip) 

Js 



g{p)\dp, 



diag{S)A j^ 

where A is the area of S, diag{S) is the diagonal of the bounding box of S. The distance is the L^ norm 
between / and g, normahzed by the diagonal of surface. In our experiment, the distance is 0.000044, 



which is very small. This shows that our quasi-conformal map method is accurate. Figure 12 shows the 
histogram of the real part, imaginary part and argument of /x/, fig and fifog. (A), (B) and (C) shows the 
histograms of the real part, imaginary part and argument of the Beltrami coefficient /x/ of / computed 
by our method. The histograms show that Real(/x/) = 0.15, Iniag(/x/) = 0.15 and arg(/x/) = 0.7854 
on almost all vertices, which agree with the exact solution. (D), (E) and (F) shows the histograms of 
the Beltrami coefficient fig of g. (G), (H) and (I) shows the histogram of the real part, imaginary part 
and argument of the Beltrami coefficient fifog of the composition map fog. The histograms show that 
Real(/x/o^) = 0.35, Imag(/x/) = 0.12 and arg(/x/) = 0.34 on almost all vertices, which agree with the 
exact solution. 

C Quasi-conformal Maps for Genus One Closed Surface 



We test our algorithm for genus one closed surface as shown in Figure [T3| We first set jjl to be zero, and 
compute a conformal flat metric using the curvature flow. Then we compute a homology basis, {a, 6} as 
shown in the leftmost frame of the top row. We then embed a finite portion of the universal covering space 
of the Kitten surface using the flat metric, shown in the right frame of the first row. The red rectangle 
shows the fundamental polygon, which is a parallelogram, with two adjacent edges Za and z^. The lattice 
r is formed by the translations generated by Za and z^, 

r = {mZa + nzi)\m^ n G Z}. 
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/ifz) = 






/4z) = -0,1-/0,1 



}a{z) =0.25 





fi{z]={)25 + t{)25 

Fig. 10 
CONFORMAL AND QUASI-CONFORMAL MAPS FOR A MULTIPLY CONNECTED DOMAIN. 



The Kitten surface can be represented as the quotient space M = ^. This gives the conformal parameter 
domain of the surface. The rightmost frame of the first row illustrates the circle packing texture mapping 
induced by the conformal parameterization. In the second and the third row, we set different Beltrami 
coefficients. /x(z) are constants for the second row. 

For the last row, the Beltrami coefficient is defined in a more complicated way. Because /x is defined 
on the Kitten surface, then it must satisfy the following consistency condition /x(z) = /x(z + mza + nz^). 
Given a point z G C, we can find a pair of real numbers a, (5 G [0, 1), such that z = aza + /3zi){modr). 
Then /i is defined as /x(z) = ^(cos27ra + zcos27r/3), which satisfies the above consistency condition. 



D. Quasi-conformal Maps for Genus Two Surfaces 

Our method can compute quasi-conformal maps for high genus surfaces, as shown in Figure [T4j We 
use hyperbolic Yamabe flow to compute the hyperbolic metric of the surface, then the homology basis 
{ai, 6i,a2, 62}- We flatten a fundamental domain, compute the Fuchsian group generators, and flatten a 
finite portion of the universal covering space of the surface. Details can be found in |l29l , ll30ll . This gives a 
conformal atlas of the surface. Because of the difference between hyperbolic metric and Euclidean metric. 



the texture mappings have seams in Figure [14] along the homology basis. Suppose z and w are two local 
parameters, differ by a Mobius transformation, then /x should satisfy the following consistency relation 
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(a) Original surface S (b) Icxhirc mapping induced by (c) (c) Ctnilormal paramcLcr (d) Texture mapping indueed by (e) 







(e) linage of / : .S — ^ D[ 
witli iJj 



(f) Texture mapping induced by (e) 



(g) Image of g : D[ — D2 
vviih Ug 



(li) Texture mapping induced by (g) 





(i) Image oi'h: S -^ D2 



(j) fexlure mapping induced by (1) 



Fig. 11 
Composed quasi-conformal maps for a topological quadrilateral. 



(See Definition 5.3): 



w 



jjiiw) — = jJiiz). 



For example, let p G S' is a point on ak p ^ ak, it has two parameters Zp G a^ and Wp e a^ , Wp = ak{zp). 
Then /x(zp) and iii{wp) should satisfy the above consistent constraint. In our experiments, we find a n-ring 
neighbor (n = 4) of ak, bk, denoted as R, then define /x(^^) = zq for all vertices Vi not in i?, fi{vi) = 0, 
for Vi in ak or 6^^. // is extended to other vertices as a complex valued harmonic function, A/x(^^) = 0, 
Vi RU ak U bk, where A is the Laplace-Beltrami operator of the original surface. This will ensure 



the consistency relation holds for /x. In Figure [14} zq is 0.2 + iO.2 for the left frame of second row, 0.3 
for the right frame of the second row, and zq = z for the last row. From the figure, we can see the 
deformation of the conformal structure of the surface (shape of the fundamental domain)with different 
Beltrami coefficients. The results show that our method can be applied effectively on general Riemann 
surfaces of high genus. 



VIII. Conclusion 

Many surface mappings are quasi-conformal in the real world. According to quasi-conformal Te- 
ichmiiller theory, in general, there exists an one-to-one map between the quasi-conformal maps and the 
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Fig. 12 
Histogram of the real part, imaginary part and argument of the Beltrami coefficients. 



TABLE I 

Computational time. 



Figure 


#vertices 


#faces 


time(s) 


M 


20184 


39984 


99 




n 




25220 


49982 


131 


1 


20184 


39984 


101 


5 


15306 


29990 


87 


10 


13515 


26304 


108 


n 


10000 


20000 


25 


n 


2057 


4118 


31 
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Fig. 13 

CONFORMAL AND QUASI-CONFORMAL MAPS FOR A GENUS ONE SURFACE. 



Beltrami coefficients. This work introduces a method to compute quasi-conformal map from Beltrami dif- 
ferentials using auxiliary metric method. The auxiliary metric is constructed from the Beltrami differential, 
such that the desired quasi-conformal map becomes a conformal one under the auxiliary metric. Theoretic 
proof shows the rigor of the method, experimental results demonstrate the generality and accuracy of 
the algorithm. Auxiliary metric method can also be applied to most existing conformal parameterization 
methods. 

APPENDIX 
PROOF OF THEOREM 4.7 



Proof: We first compute the entries in the Hessian matrix explicitly. Note that 









Since X^ = 27r — ^ 



[i,j,/c]GF 



OJ^, differentiating both sides with respect to Uj we get 



h 



dK. 



u 



39'^ del' 

duj duj duj 
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Fig. 14 

CONFORMAL AND QUASI-CONFORMAL MAPS FOR A GENUS TWO SURFACE. 



Similarly, 



i^ii 



OK, ^ 09^:' 



du 



j,k 



dui 



QQJk QQJk 

Now, -Q^ and -^^ can be computed explicitly. We will show the Euclidean case, whereas the hyperbolic 
case can be proven similarly. Firstly, we observe that by cosine rule 



2ijikCo^e, = i] + il-il 



Differentiating both sides with respect to li 



-2ljlk^in9. 



dk 



-2L 



Thus, 1^ = —li/2A, where A= area of the triangle = 2ljlksm9i Differentiating both sides with respect 



-21 j cos 6i — 2ljlk sin di 



dh 



2lu 
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Thus, g = ^^''^^^^~^' = -liCo^ej/2A, since kcosOj + IjOi = k. Also, note that g = and g: = /^ 
since L- = e^^e^^l^. We have 



Lfcj 



d9i dOi dli dOidlk , a\ioA ^a 

Similarly, 

96^i dOi dlj dOi dlk . , , ^ , , ^ \ /^ a ^ ^ 

'B~^Wb^W^^ (-/i/jCos6'fc + lilkCos6j)/2A = cotOj + cot Ok 

■ 

PROOF OF THEOREM 4.8 

Proof: The local convexity of the Euclidean Yamabe energy is due to the fact that the Hessian of 
E{u) is the matrix (hij) which is semi-positive definite. Clearly, the summation of each row is zero and 
only the diagonal elements are positive. Furthermore, since the matrix is positive definite on the linear 
space ^^Ui = 0, it follows that H is locally strictly convex on the planes. For detail, please see [36]. 

We now prove the convexity of the hyperbolic Yamabe energy. This fact was also known to Springborn 
et al. ahead of us. We prove the Hessian matrix of the hyperbolic Yamabe energy is positive definite. Let 
ai,a2,a3 be the lengths of a hyperbolic triangle. Make conformal change to produce a new hyperbolic 
triangle of lengths ?/i,?/2,2/3 so that sinh(^^/2) = sinh(a^/2)e^^+^'^, {i,j,k} = {1,2,3}. 9i represents the 
corner angle at the vertex Vi. Let H = [|^] be the matrix. 

Fact 1 : det{H) ^ for all u's. Indeed, the map from (2/1,2/2,^3) -^ {Oi, 92,0s) is a diffeomorphism. 

Fact 2 : For any (ai, a2, a3)m the set of all u = {ui, U2, us) such that yi/2 = 2sinh"'^(sinh(a^/2)e^^+^'^) 
satisfies triangular inequalities: yi + 1/2 > Vs, Vi + Vs > 1/2 and 2/2 + ?/3 > Vi form a connected set f] in M^. 

Fact 3 : Since Q is connected and H is symmetric in Q so that det{H) ^ 0, the signature of i7 is a 
constant. 

Fact 4 : Choose those ui,U2, us so that 2/1 = 2/2 = 2/3, we see easily by computation that the Hessian H 
is positive definite. 

Thus, H is positive definite over all Q. 

Now to prove Fact 2 :, introduce a new variable U = e^^+^^. The map 2/ h^ t is a diffeomorphism. 
Thus, it suffices to prove that the set Q^ = {(ti, ^2, ^3) ^ ^^iol?/^ = 2 sinh"'^(t^ sinh(^)) satisfies triangular 
inequalities } is connected. Fix ^2,^3, we will show that the set of all t such that 



I sinh '^(t2sinh( — )) — sinh '^(t3sinh( — ))| < sinh '^(t2sinh( — ))+sinh '^(t3sinh( — )) 

^ jLj Zj Zj 

is connected. This is obvious since /(t) = sinh~^ is a strictly increasing function in t. ■ 

PROOF OF THEOREM 4.9 

Proof: The solution u{f) = {u\{i), ....Ufqit)) of the discrete Yamabe flow exists for all time so 
that there are no singularities forming at time equal to infinity. This means that 2/^(t)'s are in some 
compact interval in ]R>o and also all inner angles 6^ -^ (t) are in some compact interval inside the interval 
(0,7r). The matrix (q^) = (|^) has properties that the sum of entries in every row is zero and the 
diagonal entries are negative, (qj) is symmetric and semi-negative definite. This implies that there is a 
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positive constant A so that the eigenvalues of (qj) considered as a bihnear form restricted to the subspace 

{w G ]R^|^(;l + ... + wn = 0} is always bounded by —A for all time t G [0, oc), i.e., 



J2^^j^^^j ^-^Yl 



wl 



^j 



^N 



when Y.i=i ^i = 0. 
Note that 



^ = s<2-E'^') = -E 



j:k j,k 



duj dt ^-^ duj 
3 •' j ■' 



Now, consider G{t) = J2i=ii^ii^) ~ ^if. Its derivative can be calculated as 



-xt 



We have, G\t) < -XG{t), Thus, G{t) < Ce'^^ and so 

\K,{t)-K,\<c^e-'^' 



[1 

[2 

[3 

[4: 

[5 

[6: 

[7 

[8 
[9 

[lo: 
[11 

[12 

[13 
[14; 
[15 

[16; 

[17 
[18 

[19 

[2o: 

[21 

[22 
[23 



References 

L. Ahlfors, "Conformality with respect to riemannian matrices," Ann. Acad. Sci. Fenn. Sen, vol. 206, pp. 1-22, 1955. 

H. Grotzsch, "Uber die verzerrung bei schlichten nichtkonformen abbildungen und eine damit zusammenh angende erweiterung des 

picardschen," Rec. Math., vol. 80, pp. 503-507, 1928. 

C. Morrey, "On the solutions of quasi-linear elliptic differential equations," Trans. Amer. Math. Soc, vol. 43, pp. 126-166, 1938. 

M. Lavrentjev, "Sur une classe de representations continues," Rec. Math., vol. 48, pp. 407-423, 1935. 

L. Ahlfors, Lectures in Quasiconformal Mappings. Van Nostrand Reinhold, New York, 1966. 

O. Lehto and K. Virtanen, Quasiconformal Mapping in the Plane. Springer, Berlin, 1973. 

L. Bers, "Quasiconformal mappings, with applications to differential equations, function theory and topology," Amer Math. Soc. Bull, 

vol. 83, no. 6, p. 10831100, 1977. 

L. Carleson and T. Gamelin, Complex Dynamics. Springer, New York, 1993. 

S. G. RR Belinskii and I. Yanenko, "The use of a class of quasiconformal mappings to construct difference nets in domains with 

curvilinear boundaries," USSR Comp. Math. Phys., vol. 15, pp. 133-144, 1975. 

L. Bers, Mathematical Aspects of Subcritical and Transonic Gas Dynamics. Wiley, New York, 1958. 

L. Bers and L. Nirenberg, "On a representation theorem for linear elliptic systems with discontinuous coefficients and its applications," 

Convegno Internazionale Suelle Equaziono Cremeonese, Roma, pp. 111-140, 1955. 

Bers and L. Nirenberg, "On linear and nonlinear elliptic boundary value problems in the plane," Convegno Internazionale Suelle 

Equaziono Cremeonese, Roma, pp. 141-167, 1955. 

R Daripa, "On a numerical method for quasiconformal grid generation," /. Comput. Phys., vol. 96, pp. 229-236, 1991. 

C. Mastin and J. Thompson, "Discrete quasiconformal mappings," Z Angew. Math. Phys., vol. 29, pp. 1-11, 1978. 

M. S. Floater and K. Hormann, "Surface parameterization: a tutorial and survey," in Advances in Multire solution for Geometric 

Modelling. Springer, 2005, pp. 157-186. 

A. Sheffer, E. Praun, and K. Rose, "Mesh parameterization methods and their applications," Foundations and Trends® in Computer 
Graphics and Vision, To appear. 

K. Hormann, B. Levy, and A. Sheffer, "Mesh parameterization," SIGGRAPH 2007 Course Notes 2, 2007. 

B. Levy, S. Petitjean, N. Ray, and J. Maillot, "Least squares conformal maps for automatic texture atlas generation," SIGGRAPH 2002, 
pp. 362-371, 2002. 

M. Desbrun, M. Meyer, and P. Alliez, "Intrinsic parameterizations of surface meshes," Computer Graphics Forum (Proc. Eurographics 

2002), vol. 21, no. 3, pp. 209-218, 2002. 

A. Sheffer and E. de Sturler, "Parameterization of faced surfaces for meshing using angle based flattening," Engineering with Computers, 

vol. 17, no. 3, pp. 326-337, 2001. 

A. Sheffer, B. Levy, M. Mogilnitsky, and A. Bogomyakov, "ABF-i~i-: Fast and robust angle based flattening," ACM Transactions on 

Graphics, vol. 24, no. 2, pp. 311-330, 2005. 

R. Zayer and H. P. S. B. Levy, "Linear angle based parameterization," in In Symposium on Geometry Processing, 2007, pp. 135-141. 

X. Gu and S.-T. Yau, "Global conformal parameterization," in Symposium on Geometry Processing, 2003, pp. 127-137. 



PREPRINT 24 



[24; 

[25 

[26: 

[27 

[28 
[29 
[30: 

[31 

[32 

[33 

[34: 
[35 

[36: 
[37 
[38 

[39 

[40: 

[41 

[42 

[43 

[44 
[45 
[46 
[47 



S. J. Gortler, C. Gotsman, and D. Thurston, "Discrete one-forms on meshes and applications to 3D mesh parameterization," Computer 

Aided Geometric Design, vol. 23, no. 2, pp. 83-112, 2005. 

N. Ray, W. C. Li, B. Levy, A. Sheffer, and P. Alliez, "Periodic global parameterization," ACM Transactions on Graphics, vol. 25, no. 4, 

pp. 1460-1485, 2005. 

K. P. F. Kalberer, M. Nieser, "Quadcover- surface parameterization using branched coverings." Comput. Graph., vol. 26, no. 3, p. 

375384, 2007. 

C. Gotsman, X. Gu, and A. Sheffer, "Fundamentals of spherical parameterization for 3d meshes," ACM Transactions on Graphics, 

vol. 22, no. 3, pp. 358-363, 2003. 

E. Praun and H. Hoppe, "Spherical parametrization and remeshing," ACM Transactions on Graphics, vol. 22, no. 3, pp. 340-349, 2003. 
C. Grimm and J. F. Hughes, "Parameterizing N-holed tori," in IMA Conference on the Mathematics of Surfaces, 2003, pp. 14-29. 
M. Jin, J. Kim, F. Luo, and X. Gu, "Discrete surface ricci flow." IEEE Transaction on Visualization and Computer Graphics, vol. 14, 
no. 5, p. 10301043, 2008. 

P. L. Bowers and M. Hurdal, "Planar conformal mapping of piecewise flat surfaces," In Visualization and Mathematics III (Berlin, 
2003), Springer-Verlag, pp. 3-34. 

A. I. Bobenko and B. A. Springbom, "Variational principles for circle patterns and koebe's theorem," Transactions of the American 
Mathematical Society, vol. 356, pp. 659-689, 2004. 

L. Kharevych, B. Springerborn, and P. Schroder, "Discrete conformal mappings via circle patterns," ACM Trans. Graph., vol. 25, no. 2, 
pp. 412-438, 2006. 

B. Chow and F. Luo, "Combinatorial ricci flows on surfaces," Journal Differential Geometry, vol. 63, no. 1, pp. 97-129, 2003. 

M. Ben-Chen, C. Gotsman, and G. Bunin, "Conformal flattening by curvature prescription and metric scaling," Comput. Graph. Forum, 
vol. 27, no. 2, p. 449458, 2008. 

F. Luo, "Combinatorial yamabe flow on surfaces." Commun. Contemp. Math, vol. 6, no. 5, pp. 765-780, 2004. 

A. Bobenko, B. Springborn, and U. Pinkall, "Discrete conformal equivalence and ideal hyperbolic polyhedra," In preparation. 

B. Springborn, P. Schrooder, and U. Pinkall, "Conformal equivalence of triangle meshes," ACM Transactions on Graphics, vol. 27, 
no. 3, pp. 1-11, 2008. 

C. Mastin and J. Thompson, "Quasiconformal mappings and grid generation," SI AM J. Sci. Statist. Comput., vol. 5, no. 2, pp. 305-310, 
1984. 

A. Vlasyuk, "Automatic construction of conformal and quasiconformal mapping of doubly connected and triple connected domains," 
Akad. Nauk Ukrainy Inst. Mat., preprint (Akademiya Nauk Ukrainy Institut Matematiki, preprint), vol. 57, pp. 1-57, 1991. 
J.Weisel, "Numerische ermittlung quasikonformer abbildungen mit finiten elementen," Numer. Math., vol. 35, pp. 201-222, 1980. 

P. Daripa, "A fast algorithm to solve nonhomogeneous cauchyriemann equations in the complex plane," SIAM J. Sci. Statist. Comput., 

vol. 13, no. 6, pp. 1418-1432, 1992. 

P. Daripa and D. Masha, "An efficient and novel numerical method for quasiconformal mappings of doubly connected domains," 

Numerical Algorithm, vol. 18, p. 159175, 1998. 

H. W. Guggenheimer, Differential Geometry. Dover Publications, 1977. 

H. M. Farkas and I. Kra, Riemann Surfaces. Springer, 2004. 

R. S. Hamilton, "Three manifolds with positive ricci curvature," Journal of Differential Geometry., vol. 17, pp. 255-306, 1982. 

B. Chow, "The ricci flow on the 2-sphere," /. Differential Geom., vol. 33, no. 2, pp. 325-334, 1991. 



